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f**^ ' Critical velocities of supcrfluid Fermi gases in optical lattices are theoretically investigated across 

, the BCS-BEC crossover. We calculate the excitation spectra in the presence of a supcrfluid 

^■f^ , flow in one- and two-dimensional optical lattices. It is found that the spectrum of low-lying 

' Anderson-Bogoliubov (AB) mode exhibits a roton-like structure in the short-wavelength region 

• ' due to the strong charge density wave fluctuations, and with increasing the superfluid velocity 

one of the roton-like minima reaches zero before the single-particle spectrum does. This means 
that superfluid Fermi gases in optical lattices are destabilized due to spontaneous emission of the 
roton-like AB mode instead of due to Cooper pair breaking. 
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H 1. Introduction a new question, i.e., how does the criti- 

. J. .Ti -i- a -J cal velocity for bosonic superfiuids change 

According to Landau criterion, supernuid •' f b 

n r u • 4-- 1 i. -ii, in the BCS regime? Recently, Miller et al. 

now oi bosonic particles propagates with- 

1-4- 1 -i investigated critical velocities of superfluid 

out dissipation, unless its group velocity ex- ° 

J i • -i- 1 1 J i Fermi gases across the BCS-BEC crossover 

ceeds a certain critical value and spontaneous ° 

.. r ■ ■ J ji rru J- in a moving one-dimensional (ID) optical 

creation oi excitations is induced. Ihe di- ° \ j f 

. . , , r ii, lattice. They determined the critical veloci- 

rect experimental verincations of the crite- 

, , , , , J • 1 • ties at which the abrupt decrease of fermionic 

rion have been attempted mainly m super- ^ 

fl-j4xj2 jj-w + ■ T3 3 pairs in the condensate occurs, when the ve- 

nuid He and dilute atomic Bose gases. ^ ' 

rpi + J- I, ■ ■ A ■ A- locity of the lattice potential is varied, 

ihese studies have inspired various discus- ^ 

A A ^ A- ( Theoretically, critical velocities in the 

sions and contributed to understanding oi 

;i a -A \, ^f BCS-BEC crossover have been investigated 
the supernuid phenomenon itseli. 

The recent realization of the BCS-BEC ^ ^^^^^^ ^y^*^"^' '^l™^ ^'^^^ 

crossover in ultracold Fermi gases^ has raised instability of superfluid flow in a uniform 
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system is induced by Cooper pair breaking in 
the BCS regime.® In this paper, motivated 
by the experiment done by Miller et al.,^ 
we investigate critical velocities of superfluid 
Fermi gases in the BCS-BEC crossover in op- 
tical lattices. We calculate excitation spectra 
and determine critical velocities in ID and 
2D optical lattices. We find that the exci- 
tation spectrum of the low-lying Anderson- 
Bogoliubov (AB) modc®'^'^ has a character- 
istic roton-like structure, and the instability 
of superfluid flow in the BCS regime is in- 
duced by a different mechanism from that in 
a uniform system: the spontaneous emission 
of AB mode destabilizes the superfluidity. 

2. Formalism 

We consider two-component atomic Fermi 

gases with equal populations. We as- 
sume that the lattice potential is sufficiently 
deep that the tight-binding approximation is 
valid. Thus, the system can be described by 
a single-band Hubbard model 

H = -J (4cja + /i.e.) 

+uY^cl^c\^CiiCi^, (1) 

i 

where Cj^ is an annihilation operator of an 
atom with mass m and pseudospin a =1, J, on 
the j-th site, J is a hopping energy between 
nearest-neighbor sites, and U{< 0) is an on- 
site attractive interaction. We set h = 1 
throughout the paper. 

We assume that a superfluid flows with 
the velocity v in the coordinate system fixed 
with respect to the lattice potential, i.e. 
Cooper pairs have the center of mass quasi- 
momentum q = 2mv. The corresponding 
order parameter is given by 

= —^'Y^{C'-k+mviCk+mv'[)^ (2) 
k 

where M is the number of lattice sites. 
By diagonalizing the Hamiltonian within the 



BCS mean-field approximation, we obtain 
the single-particle excitation energy as 

Et{k) = 77„(fe) ± V(£t,(fe)-M)2 + |A,|2, 

(3) 

where rjv{k) — 2JTi^ sm.{mv^d) sin^k^d) and 
e„(fe) = 2JEi,(l — cos{mVud) cos{kvd)) with 
d being the lattice constant. The en- 
ergy gap and the chemical potential ii 
are determined by solving the gap equation 
and the number equation self-consistently. 
This scheme interpolates the BCS and BEC 
regimes at T = 0.^ 

To determine the critical velocities, it is 
necessary to calculate the excitation spectra 
in the presence of a superfluid flow. For neu- 
tral fermionic superfluids, there exists a low- 
lying AB mode®''"'® in addition to the single- 
particle excitation. The spectrum of the AB 
mode can be obtained from the pole of the 
density response function x(9, w).^'^ We cal- 
culate x(9)'^) by employing the generalized 
random phase approximation (GRPA) devel- 
oped by Cote and Griffin.^" 

3. Results 

Conditions for Cooper pair breaking are de- 
termined from the single-particle excitation 
energy Eq. (3). When the absolute value of 
the superfluid velocity v exceeds a certain 
critical value Upb, one flnds that E:^{k) be- 
comes negative. The negative energy of the 
single-particle excitation indicates the occur- 
rence of spontaneous Cooper pair breaking. 
Thus, the critical velocity for Cooper pair 
breaking is given by Vp^.^'^'^^ 

In order to compare the instabilities in- 
duced by two kinds of excitations, i.e., the 
single-particle excitation and AB mode ex- 
citation, it is necessary to monitor both the 
excitation spectra in parallel. For that pur- 
pose, it is useful to plot the dynamical struc- 
ture foctor S{q,Lv) = —lmx{q,io)/n. First, 
we investigate the case of ID optical lattices. 
In Fig. 1, S{q,u)) is shown for (a) current- 
free and (b) current-carrying cases. In Fig. 1, 
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Fig. 1. Dynamical structure factors S{q,u>) in ID 
optical lattices are shown for (a) current-free {v = 0) 
and (b) current-carrying cases {v = 0.095/md), when 
n = 0.25 and U = — l.OJ. The dark region indicates 
large S{q,u!). In each panel, the upper gray region 
corresponds to the particle-hole continuum, and the 
lower dashed curve is the spectrum of the AB mode. 



the upper gray regions are the particle-hole 
continuums, and the lower boundaries of the 
particle-hole continuums correspond to u = 
min [E:^{k + q) — E~{k)]. The lower curves 
in Fig. 1 are the spectrum of AB mode, 
which becomes phonon-like when qd ^ 1.^ 
It should be noted that the spectrum of the 
AB mode has roton-like minima at \q\ ~ 2fci? 
where kp is the Fermi wavevector. Such 
a roton-like structure also appears in the 
spectrum of AB mode in uniform ID Fermi 
gases. The energies of the roton-like min- 
ima decrease as one approaches the half- 
filling when w = 0. They eventually become 




Fig. 2. Critical velocities Vp^, and ii^b as functions 
of the filling n when U = —l.OJ. The solid line 
represents r>pb and the dashed line represents i^ab- 



zero at fc = Tr/d at the half- filling. This sug- 
gests that exciting density fluctuations with 
k = n/d costs zero energy and the super- 
fluid ground state is in competition with the 
charge density wave (CDW) state. Thus, the 
roton-like minima reflect the strong CDW 
fluctuations. 

In the presence of a superfluid flow with 
V > 0, the spectrum slants towards the left 
side (see Fig. 1(b)). The critical velocity due 
to pair breaking Wpb coincides with the ve- 
locity at which the particle-hole continuum 
touches the line ui — 0. Since the AB mode 
spectrum locates below the particle-hole con- 
tinuum, one of the roton-like minima reaches 
zero at the velocity Wab which is smaller than 
Wpb- According to the Landau criterion,^ this 
means that the spontaneous emission of the 
AB mode is induced before the Copper pair 
breaking occurs. In Fig. 2, the critical ve- 
locities Upb and Wab are plotted as functions 
of the filling n. One clearly sees that Vab 
is always smaller than Vpb, and that both 
Wab and Vph decrease monotonically with in- 
creasing filling. The difference between the 
two velocities diminishes near the half-filling 
(n^l). 

We proceed to discuss Fermi gases in 2D 
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Fig. 3. Dynamical structure factor S{q,Lo) in 2D 
optical lattices for qx = qy. We set U = —4. 5 J, 
n = 0.5, and v = {v, v)/^j2 with v = 0A67/md. The 
dark region indicates large S{q,ijj). The upper gray 
region corresponds to the particle-hole continuum, 
and the lower curve is the spectrum of AB mode. 

optical lattices. In Fig. 3, we plot the dy- 
namical structure factor for = Qy, 
when V = (u/\/2, v/^/2). Similarly to the ID 
case, the spectrum of AB mode has a roton- 
like minimum in Fig. 3. It reaches zero be- 
fore the particle-hole continuum does when 
V is increased. Thus, Uab is smaller than iipb, 
and the instability induced by emission of the 
AB mode also dominates in 2D optical lat- 
tices. Furthermore, it is found that the di- 
rection of the superfluid velocity v affects the 
spectrum of AB mode. The energy of the 
roton-like minimum decreases more rapidly 
for V = {v/V2,v/V2) than for v = {v,0). 
This implies that the CDW fluctuation is en- 
hanced when the supercurrent is flowing in 
the direction of (tt, tt). 

4. Conclusion 

In summary, we have studied critical veloci- 
ties of superfluid Fermi gases in ID and 2D 
optical lattices using the attractive Hubbard 
model. We calculated the excitation spectra 
of single-particle excitation and AB mode by 
employing the GRPA. We have shown that 
the spectrum of AB mode has a roton-like 



structure at \q\ ~ 2fci?, and it always lies be- 
low the particle-hole continuum due to the 
strong CDW fluctuations in optical lattices. 
As the velocity of the superfluid flow is in- 
creased, the energy of a roton-like minimum 
decreases and the instability due to the emis- 
sion of roton-like excitations was shown to 
occur before Cooper pair breaking occurs. 
We note that Burkov and Paramekanti have 
proposed another type of current-induced in- 
stability of superfluid Fermi gases in optical 
lattices. 

Acknowledgments 

We would like to thank K. Kamide, S. Kuri- 
hara, and T. Nikuni for valuable comments 
and discussions. I.D. and D.Y. are supported 
by a Grant-in- Aid from JSPS. 

References 

1. L. D. Landau, J. Phys. (USSR) 5, 71 (1941). 

2. D. R. Tilley and J. Tilley, Superfluidity and 
Superconductivity, (Bristol, Hilger 1991). 

3. R. Onofrio et al, Phys. Rev. Lett. 85, 2228 
(2000). 

4. C. A. Regal et al, Phys. Rev. Lett. 92, 
040403 (2004); M. W. Zwierlein et al, ibid. 
92, 120403 (2004). 

5. D. E. Miller et al, Phys. Rev. Lett. 99, 
070402 (2007). 

6. R. Combescot et al, Phys. Rev. A 74, 042717 
(2006). 

7. P. W. Anderson, Phys. Rev. 112, 1900 
(1958). 

8. L. P. Pitaevskii et al, Phys. Rev. A 71, 
053602 (2006). 

9. A. J. Leggett in Modern Trends in the Theory 
of Condensed Matter, edited by A. Pekalski 
and R. Przystawa (Springer- Verlag, Berlin, 
1980). 

10. R. Cote and A. Griffin, Phys. Rev. B 48 
10404 (1993). 

11. M. Rodriguez and P. Torma, Phys. Rev. A 
69, 041602(R) (2004). 

12. T. Aim and P. Schuck, Phys. Rev. B 54, 
2471 (1996). 

13. A. A. Burkov and A. Paramekanti, Phys. 
Rev. Lett. 100, 255301 (2008). 



